Quantum criticality and universal scaling of strongly attractive spin-imbalanced Fermi 

gases in a ID harmonic trap 
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We investigate quantum criticality and universal scaling of strongly attractive Fermi gases con- 
fined in a one-dimensional (ID) harmonic trap. We demonstrate from the power-law scaling of the 
thermodynamic properties that current experiments on this system are capable of measuring univer- 
sal features at quantum criticality, such as universal scaling and Tomonaga-Luttinger liquid physics. 
The results also provide insights on recent measurements of key features of the phase diagram of 
a spin-imbalanced atomic Fermi gas [Y. Liao et al, Nature 467, 567 (2010)] and point to further 
study of quantum critical phenomena in ultracold atomic Fermi gases. 
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Very recently, the one-dimensional (ID) strongly at- 
tractive two-component Fermi gas has attracted much at- 
tention from theorists [lH7| and experimentalists [8| due 
to the existence of an exotic pairing mechanism. Inves- 
tigation [ij-HH shows that this novel pairing is closely 
related to the elusive Fulde-Ferrel-Larkin-Ovchinnikov 
(FFLO) [HI state involving BCS pairs with nonzero 
center-of-mass momenta. 

The ID Fermi gas is one of the most important exactly 
solvable quantum many-body systems. It was solved long 
ago by Yang [l3[ and Gaudin 14 1 using the Bethe ansatz. 
Although the study of the attractive Fermi gas was initi- 
ated soon after [15l . [16| , it was not until much later that 




this model began to receive more attention [17[ • In terms 
of the polarization p the model exhibits three quantum 
phases at zero temperature the fully paired phase 

which is a quasi-condensate with P = 0, the fully polar- 
ized normal Fermi liquid with P = 1, and the partially 
polarized (ID FFLO-like) phase for < P < 1. 

For a trapped imbalanced Fermi gas it is found [l|, 0] 
within the local density approximation (LDA) that a par- 
tially polarized ID FFLO-like state sits in the trapping 
center surrounded by wings composed of either a fully 
paired state or a fully polarized Fermi gas. The phase 
boundaries of the zero-temperature phase diagram shown 
in Fig. [1] can be determined precisely from the exact solu- 
tion by the vanishing of the axial density difference (solid 
line) and the minority state axial density (dashed line). 
The key features of this T = phase diagram were ex- 
perimentally confirmed using finite temperature density 
profiles of trapped fermionic 6 Li atoms [8|. 

Most recently, quantum criticality and universal scal- 
ing behaviour 18[ are being explored in low-dimensional 
cold atomic matter. For example, in experiments on the 
2D Bose gas [l9| , following a theoretical scheme for map- 
ping out quantum criticality 0, 21 1- From this view- 



FIG. 1: Phase diagram for two-component fermions in an har- 
monic trap at zero temperature as a function of polarization 
with Na\ D /a? z — 0.2. Here FP denotes the fully paired phase, 
F the fully unpaired phase and PP the FFLO-like phase. V 
is the vacuum phase. 



phase transitions at T — 0, is particularly valuable to 
test universal scaling through finite temperature density 
profiles of trapped ultracold atoms (22|. Here we illus- 
trate that finite-temperature properties of the quasi-lD 
trapped Fermi gas allow the exploration of a wide range 
of physical phenomena, such as universal Tomonaga- 
Luttinger liquid (TLL) physics, scaling theory and the 
nature of the FFLO state at quantum criticality. 

The Model and equation of state.- We consider the ID 
^-interacting attractive spin-1/2 Fermi gas with N — 

+ fermions of mass m and external magnetic field 
H. The system is described by the hamiltonian 
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in which the three terms are kinetic energy, interaction 
energy and Zeeman energy, respectively. Here the inter- 
component interaction is related to an effective ID scat- 
tering length Qi n = — — — which can be tuned from 
the weakly interacting regime (gm — > - ) to the strong 
coupling regime (gin — oo) via Feshbach resonances 
and optical confinement. For convenience, we define the 
interaction strength as c = mguj/h 2 and dimensionless 
parameter 7 = c/n for physical analysis, where n = N/L 
is the linear density and L is the length for system. We 
set Boltzmann constant fes = 1. 

For the strongly attractive spin-1/2 Fermi gas at finite 
temperatures, the thermodynamics of the homogeneous 
system is described by two coupled Fermi gases of bound 
pairs and excess fermions in the charge sector and ferro- 
magnetic spin-spin interaction in the spin sector. Spin 
fluctuations are suppressed by a strong effective mag- 
netic field at low temperatures. For the physically in- 
teresting low temperature and strong coupling regime, 
i.e., T <C Cb, H and 7 » 1 a high precision equation 
of state [22h24| can be derived from the thermodynamic 
Bethe ansatz (TBA) equations [25| in terms of the Yang- 
Yang grand canonical ensemble [26(. Using the binding 
energy Eb = h 2 c 2 /Am as the unit of energy and defining 
/2 = /i/ffh, h = H /eb, t = T/eb, p = p/ \ceb\, the pressure 
p = p b + p u of the system is found to be 
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where /* = Li„ (-e A ^ 1 ) and /» - Li* 
of the standard polylog function Li„ (x), with 
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fc= 



f s = te- h / t e- 2 P"/ t I (2p u /t) and I n (x) = 
k\{n+k)\ {x /2) n+2k . The thermodynamic quantities 
such as the particle density n, n u for unpaired fermions, 
n b for paired fermions, and compressibility follow from 
Eq. ([2]) and the standard thermodynamic relations. The 
total pressure serves as the equation of state which pro- 
vides high precision thermodynamics over a wide tem- 
perature range T < 0.2Eb, see (22I. l23j. 

Quantum phase diagram in a harmonic trap. - For spin 
imbalanced attractive fermions in an harmonic trap, the 
equation of state © can be reformulated within the LDA 
by the replacement [i (z) = \i (0) — \mui 2 z z 2 in which z 



is the position and u> z is the frequency within the trap. 
Using the dimensionless chemical potential this becomes 



jj,(z) = /t(0) — 2z 2 , where z = z/ 
monic characteristic radius a z 
tal particle number N = J_ n (z) dz and polarization 
P = n u (z) dz/N become 



\c\) with the har- 
= \Jfij (mu} z ). The to- 



Na\ D /a 2 = 4 / h(z)dz 
J —00 

/oo 
h u (z)dz x a 2 z l (Na 2 1D ) (4) 
-00 

where h(z) = 1/ \j (z)\ and n u (z) = n u (z) /|c|. 

At finite temperatures, the phase boundaries can be 
determined from the equation of state © within the LDA 
((3]) , where the boundaries of vanishing density difference 
(black solid line) and vanishing unpaired fermions (red 
dashed line) form three phases (see Fig.[T]). As t — > 0, the 
phase boundaries determined via ((2]) are consistent with 
the zero-temperature results. For strong attraction with 
polarization, the atoms with opposite spin states form 
hard-core bosons which are strongly correlated with ex- 
cess fermions. The polarization can be changed by tuning 
the effective magnetic field. This results in two extreme 
phases - the fully paired and fully polarized phases. At 
zero temperature, these phases intersect at global polar- 
ization P = P r , where we find 
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Here the phase boundaries are sharp lines and P c is the 
intersection point in Fig. [TJ At finite temperatures the 
intersection point P c decreases as the temperature in- 
creases. This result is consistent with experimental ob- 
servation However, for the experimental temperature 
t = 175 nK (about 0.0267£&), our theoretical intersection 
point is smaller than the experimentally estimated value. 
This discrepancy originates merely in defining the edge 
of the cloud in the experiment where the phase bound- 
aries are no longer sharp lines 27J. We also find that in- 
creasing temperature results in shrinking the fully-paired 
phase and populating the vacuum. The existence of these 
quantum phases at low temperatures is a manifestation 
of different TLL phases, i.e., a two-component TLL of 
FLLO-like states and a TLL of hard-core bosons and a 
TLL of unpaired fermions (23j , where the low-lying exci- 
tations are close to the Fermi points. 

The FFLO-like state behaves like a mixture of com- 
posite fermions with mass 2m and excess fermions with 
mass m [22j]. The density of state changes critically as 
the driving parameter chemical potential is varied across 
the phase boundaries. Fig. [2] shows the density distri- 
butions for Na\ D /a 2 — 0.2 at zero temperature and for 
two different values of polarization. For small polariza- 
tion (10%), an FFLO-like state coexists in the trap center 
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FIG. 2: The zero-temperature density distribution profiles for 
trapped fermions with Na\ D /a 2 z = 0.2 for (a) P = 10% and 
(b) P = 50%. Here we n u — n\ — n± and n b — nj.. 



with a BCS-like state at the edges; while for large polar- 
ization (50%), the FFLO-like state lies in the trap center 
with a fully polarized normal state at the edges. We 
also find the FFLO-like state occupies the whole trap for 
the critical polarization (16%), which is consistent with 
the experimental observation [8]. We now demonstrate 
how to map out the zero-temperature phase diagram and 
quantum criticality from finite temperature density and 
compressibility profiles of the trapped gas. 
Quantum criticality in a harmonic trap.- The system ex- 
hibits universal scaling behaviour in the vicinity of the 
critical points as a result of the many-body effects. The 
thermodynamic functions of the homo gene ous gas can be 
cast into a universal scaling form Tj^, 2^, 21 1. E.g., the 
density and compressibility scale as 



n(»,T) = nQ + T i^g(^-y (6) 

«( M ,T) = KQ + T^-^T^^y (7) 

Here the dimensionality d = 1, dynamical critical expo- 
nent z = 2 and correlation length exponent v = 1/2 for 
the ID strongly attractive Fermi gas 22]. The universal 
scaling functions Q(x) and IF(x) describe quantum and 
thermal fluctuations at quantum criticality after subtrac- 
tion of the background density no and compressibility kq . 

For small polarization P < P c , the chemical potential 
passes the lower critical point n C 2 — — \ from the vacuum 
into the fully paired phase then passes the upper critical 
point u c4 fa -§ + ^(1 - h)i + ^(l - h) 2 from the fully 
paired phase into the FFLO-like phase. At finite tem- 
peratures, contour plots of the entropy clearly indicate 
universal critical behavior near the critical points, see 
Fig. [3Ja). The typical T^-shape crossover temperatures 




FIG. 3: Quantum criticality for low polarization, (a) contour 
plot of entropy in the t — fi plane where t = T/sb- Dia- 
monds and triangles indicate crossover temperatures separat- 
ing the quantum critical regimes from vacuum, single compo- 
nent TLLp of paired fermions and two-component TLLpp of 
FLLO-like states. The density profiles demonstrate how to 
locate the critical points (b) /i C 2 and (c) ^t C 4 and thus map 
out their phase boundaries. The corresponding compressibil- 
ity curves (d) and (e) intersect at the critical points after a 
subtraction of the background compressibility. 




FIG. 4: Quantum criticality for high polarization, (a) con- 
tour plot of entropy in the t — fi plane where t — T/et- Dia- 
monds and triangles indicate crossover temperatures separat- 
ing the quantum critical regimes from vacuum, single compo- 
nent TLL f of unpairs and two-component TLL pp of FLLO- 
like states. The density profiles demonstrate how to locate 
the critical points (b) /i c i and (c) fj, C 3 and thus map out their 
phase boundaries. The corresponding compressibility curves 
(d) and (e) intersect at the critical points after a subtraction 
of the background compressibility. 
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separate the quantum critical regimes of the hard-core 
bosonic TLLp phase and the two-component FLLO-like 
TLLpp phase near the critical points \x C 2 and p, c 4, re- 
spectively. Except for the left-most line separating the 
vacuum from the critical regime, all of the lines can be 
determined by the breakdown of the TLL, i.e., when the 
entropy curves deviate from linear T-dependent relations 
for fixed values of /x and h. These boundaries indicate a 
crossover from linear dispersion into nonrelativistic dis- 
persion rather than phase transitions at finite T. 

For large polarization P > P c , the chemical potential 
passes the lower critical point /u c i = —h/2 from the vac- 
uum into the fully unpaired phase then passes the upper 
critical point ^ c3 « -\ (l - -^(h - l)i - ^{h - x ) 2 ) 
from the fully unpaired phase into the FFLO-like phase, 
see Fig. |U[a). Here the typical V -shape crossover tem- 
peratures indicate quantum criticality of the largely po- 
larized Fermi gas near the two critical points \i c \ and /x C 3, 
respectively. 

The zero temperature phase boundary separating the 
vacuum from the fully paired phases can be mapped out 
from the universal scaling behaviour of the pair den- 
sity curves for different temperatures, which intersect 
at the critical point /i C 2, see Fig. [3]Jb). Similarly the 
phase boundary separating the fully paired phase from 
the FFLO-like phase can be mapped out from the den- 
sity profiles of unpaired fermions in the trapped gas at 
different temperatures, i.e., the unpaired density curves 
for different temperatures intersect at the critical point 
/i C 4, see Fig. [3jc). Here we have chosen the den- 
sities without background near the critical points for 
practical purposes. We see that the density curves for 
temperatures below 0.026;, well intersect at the critical 
point. Accordingly, the compressibilities at low temper- 
atures also well intersect at the critical points, see Fig. 
(3]Jd,e). Here the background compressibility near [i c ^ is 
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The results demon- 



strate that universal sca ling behaviour persists in the 
trapped cloud of atoms [28[. As the temperature de- 
creases, we see that the compressibility evolves a round 
peak sitting in the phase of the higher density of state. 
It diverges at zero temperature. 

For the high polarization case, the density profiles of 
unpaired and paired atoms can be used to map out the 
phase boundaries \x c \ (V — > F ) and fi C 3 {F — > PP), re- 
spectively, see Fig. |UJb,c). Moreover, we show that the 
corresponding compressibility curves of the trapped gas 
at low temperatures intersect at the critical points yU c i 
and fics, see Fig. [3jd,e) . Near ji c ^ the background com- 



pressibility is k of ^ = - _ x . 

the universal scaling for other thermodynamic quantities 
such as the specific heat and magnetization are testable 
through the quantum criticality of the trapped gas. This 
provides a reliable way to determine quantum phase dia- 
grams and test universal scaling theory and TLL physics 



in ID quantum gases of cold atoms. 
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